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ABSTRACT 

The spatial distribution of Lya forest absorption systems towards ten QSOs has been 
analysed to search for large-scale structure over the redshift range 2.2 < z < 3.4. The 
QSOs form a closely spaced group on the sky and are concentrated within a 1 deg 2 
field. We have employed a technique based on the first and second moments of the 
transmission probability density function which is capable of identifying and assessing 
the significance of regions of over- or underdense Lya absorption. We find evidence for 
large-scale structure in the distribution of Lya forest absorption at the > 99 per cent 
confidence level. In individual spectra we find overdense Lya absorption on scales of 
up to 1200 km s" 1 . There is also strong evidence for correlated absorption across line 
of sight pairs separated by < 3 h" 1 proper Mpc (qo = 0.5). For larger separations the 
cross-correlation signal becomes progressively less significant. 

Key words: large-scale structure of Universe - intergalactic medium - quasars: 
absorption lines 
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1 INTRODUCTION 



([Lanzetta et al. 1995; |Le Brun, Bergeron, fc Boisse 1996 



The Lyp fuiesf seen hi Llie specQd uf disLaiiL QSOs may luii- 
oubatantial fraction of the baryonie content of the 



stitute 

Universe (Rauch & Haehnelt 1995) and its evolution can 
be traced over most of the history of the Universe. Every 
QSO absorption spectrum provides us with a representa- 
tive, albeit one-dimensional sample of the baryonie matter 
distribution, unaffected by any luminosity bias. The inti- 
mate relationship of the absorbing gas with the thermal, 
chemical, and dynamical history of the Universe makes the 
Lya forest a versatile, but not always user-friendly tool in 
the study of cosmology. 



Thb observational and theoretical advances of the pasl 
decade ^rave significantly altered our understanding of the 
Lya forest. The advent of HST's UV spectroscopic capa- 
bilities led t o the first detailed analyses of the low r edshift 
Lya forest ([Morris et al. 199l|; |Bahcall et al. 199l|). Sub- 



sequently, using multislit spectroscopy and deep imaging ot 
galaxies in the fields of QSOs, significant numbers of coinci- 
dences between the redshifts of absorption lines and galax- 
ies were found at z ^ 1. Several groups determined that 
galaxies have absorption cross-sections of ~ 200 h" 1 kpc 
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Bowcn, Blades, fe Pcttini 1996ft . Moreover, Lanzetta et al 



(1995| ), |Chen et al. (19980 , and flripp, Lu, fc Savage (1998 



found an anti-correlation of the Lya rest equivalent width, 
W, with the distance of the absorbing galaxy to the line of 
sigh t of the QSO as we ll as a s light correlation wi th luminos- 
ity (|Chen et al. 1998J). Thus |Chen et al. (1998J) confirmed 



earlier results that the extended gaseous halos of galaxies 
are directly responsible for a significant fraction of the Lya 
forest at z is 1. However, there is evidence that this may 
only be true for stronger lines (W <; 0.3 A) and that the 
weaker absorption lines trace the large-scale gaseo us struc- 
tures in which ga l axies are presumably embe dded (Le Brun 
fc Bergeron 1998J; [Tripp, Lu, fe Savage 199§). 



At high redshift, the discovery of measurable amounts 
of C IV associated with 75 per cent of all L ya absorbers 
with coloumn density A f (H I) > 10 14 ' 5 cm" 2 (Cowie et al 



1995; Songaila & Cowie 1996) in high-quality spectra ob- 
tained with the HIRES spectrograph on the Keck 10-m tele- 
scope has challenged the original notion of the Lya fores t 
arising in pristine, primordial gas. However, Lu et al. (1995) 



found almost no associated C IV absorption for systems with 
N(U I) < 10 14 cm" 2 and derived [ C/Hl < -3.5, about a fac - 



tor of te n smaller than infe rred by Songaila & Cowie (199£) 
(see also Dave et al. 1998a), suggesting a sharp drop in the 
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metalli city of the Lya forest at N(H I) ~ 10 cm . Never- Chernomordik 1995; Cristiani et al. 1995; Kulkarni et al 



theless, the metallicity of the higher column density systems 



199C ; Cristiani et al. 1997; Khare et al. 1997), and some in- 



raises the possibility of an association of these systems with 
galaxies at high redshift. Using the same C IV lines to re- 
solve the structure of the co rresponding blended Lya lines, 
Fernandez-Soto et al. (1996) have shown that the clustering 



properties of these systems are consistent with the clustering 
of present-day galaxies if the correlation length of galaxies 
is allowed to evolve rapidly with redshift (e « 2.4). 

On the theoretical side, models have progre s sed from 



vestigators found considera bly stronger signals (Ulmer 199C ; 
Fernandez-Soto et al. 1996). There have also been detections 
on larger scales. Meiksin & Bouchet (1996) used a nearest 
neighbour st atistic to derive correl ations on scales of 0.5- 
3 h _1 Mpc. Pando & Fang (1996) employed the discrete 



the early pressure-confined (Sargent et al. 198C; Ostrikei 



& Ikeu chi 1983J |lkeuchi fc Ostriker 1986} [Willigcr fc Babul 



1992) and dark matter mini-halo (|Rees 1986 Ikeuchi 198£) 



scenarios to placing the Lya forest fully within the context 
of the theory of CDM dominated, hierarchical structure for- 



mat ion. Both semi-analytic al ( Petitjean, Mticket, fc Kate: 



199E|; Bi fc Davidsen 1997} |Hui, Gnedin, fc Zhang 1997} 



Gnedin fc Hui 1998[) and full hydrodynamical numerical sim- 



ulations (|Ccn ct al. 1994J; Zhang, Anninos, fc Norman 1995; 



Miralda - Escude et al. 1996; pernquist et al. 1996 ; Wadslcy 
& Bone 1997} Theuns et al. 1998[ ) of cosmological structure 
formation, which include the effects of gravity, photoionisa- 
tion, gas dynamics, and radiative cooling, have shown the 
Lya forest to arise as a natural by-product in the fluctuating 
but continuous medium which forms by gravitational growth 
from initial density perturbations. The simulations are able 
to match many of the observed proper ties of the Lya absorp- 



wavelet transform to demonstrate the existence of clusters 
on scales of 10-20 h _1 Mpc in the Lya forest at a significance 
level of 2-4<r and that the number density of these clusters 
decreases with increasing redshift. More recently, they ob- 
served non-Gaussian behaviour of Lya forest lines on scales 
of 5-10 h -1 Mpc at a c onfidence level l arger than 95 per cent 
( |Pando fc Fang 1998; ). |Mo et al. (1992b even reported 60 and 
130 h~ Mpc as characteristic scales of the Lya forest. All 
of these studies however are based on analyses of individual 
lines of sight. 

By comparing the absorption characteristics in two (or 
more) distinct, but spatially close lines of sight, it is possi- 
ble to investigate directly the real space properties of Lya 
systems on various scales. The multiple images of gravita- 
tionally lensed QSOs have been used to establish firm lower 
limits on the sizes of Lya clouds. Smette et al. (1995) found 



uns, Leonard, fc Efstathiou 1998). A common feature of all 



tion to within reasonable accuracy ( Miralda-Escude et al. 
1996:; fchang et al 1997J [Dave et al. 1997|; |Dave et al. 199811 



Miicket et al. 1996|; |Ricdigcr, Petitjean, fc Miickct 1998| ; |The- 



that the Lya forests of the two images of HE 1104-1805 
(separation 3.0 arcsec) are virtually identical and inferred a 
2<r lower limit of 50 h _1 kpc on the diameter of Lya clouds, 
assumed to be spherical, at z — 2.3. More information on the 
sizes of Lya absorbers has been gained from the studies of 
close QSO pairs with separations of several arcsec to ~ 7 ar- 
cmin. All of the most recent analyses have concluded that 
Lya absorbers have diameters of a few hundred kpc (Fane 



the sim rlations is that the absorbing structures exhibit a va- et al. 1996 Dinshaw ct al. 1997 ; Dinshaw et al. 1998 ; Petit 

riety of 



geometries, ranging from low density, sheet-like and 



filamen ;ary structures to the more dense and more spherical 



regions where the filaments interconnect and where, presum- 
ably, galaxies form. The dividing line between these different 
geometries lies in th e range 10 14 cm" 2 Jj N 5s 10 15 cm -2 



(Cen fc Simcoe 1997; Zhang et al. 199£ 



Thus a consistent picture may be emerging: absorption 
lines with iV <; 10 14 ' 5 cm -2 are closely associated with the 
large, spherical outer regions of galaxies, clu ster strongly on 
small v eloc ity scales along the line of sight ( Fernandez-Sotc 
et al. 1996), and have been contaminated with metals by 



jean et al. 1998] ; |D'Odorico et al. 1998J ) and |Crotts fc Fang 
(1998| ) found that correlations among neighbouring lines of 
sight persist for separations up to 0.5-0.8 h -1 Mpc for lines 
with W > 0.4 A. A tentative detection of increasing cloud 
size with decreasing redshift (as would be expected for ab- 
sorbers expanding with the Hubble flow) was reported by 
Fang et al. (1996| ) and [Dinshaw et al. (1998[ ). The possi- 
bility remains however, that the coincidences of absorption 
lines are due to spatial clustering of absorbers and that the 
inferred 'sizes' in fact are an indication of their correlation 



supern ovae from a po stulated Population III (e.g. Miralda- 
Escude & Rees 1997) or by galaxy mergers (e.g. Gnedin 



length (Dinshaw et al. 1998; Cen & Simcoe 1997) as may be 



199£ ) ; ^ /hereas lower column density lines trace the intercon- 
necting, filamentary structures of the intergalactic medium. 

Whatever the case may be, it seems likely that the large- 
scale distribution of the Lya absorption holds important 
clues to its origin. Lya absorbers are probably fair tracers 
of the large-scale cosmic density field and should thus be 
able to constrain structure formation models. To investigate 
further the connection of the Lya forest with cosmological 
structure we study in this paper its large-scale structure 
both in velocity and real space by considering a close group 
of ten QSOs. The QSOs are contained within a ~ 1 deg 2 
field so that the Lya forest is probed on Mpc scales. 

There have been ma ny studies of the clustering prop- 
erties of the Lya forest. Webb (1987 ) was the first to re- 
port a weak signal in the two-point correlation function of 
fitted absorption lines on scales of ^ 100 km s -1 . This re- 



evidenced by the correlation of the estimated ' size' with line 
of sight separation found by Fang et al. (1996J ) . 

Analyses to determine the shape of absorbers , as pro- 
posed e.g. by Charlton, Churchill, & Linder (1995), have so 
far been inconclusive. Several authors agree that the current 
data are incompatible with uniform-sized spherical clouds 
but are unable to decisively distinguish between spherical 
clouds with a distr ibution of sizes, fl attened disks, or fil- 
aments and sheets (Fang et al. 1996; Dinshaw et al. 19971; 



suit was later confirmed (Ostriker, Bajtlik, fc Duncan 1988 



D'Odorico et al. 1998). Howe ver, from the analysis of a QSO 
triplet |Crotts fc Fang (199S ) found evidence that lines with 
W > 0.4 A arise in sheets. 

For still larger line of sight sep arations most work ha s 
concentrated on metal absorption. Williger et al. (1996), 
e.g., used a group of 25 QSOs contained within a ~ 1 deg 2 
region to identify structure in C IV absorption on the scale 
of 15-35 h _1 Mpc. Within this group, a subset of ten QSOs 
is suitable for studying the large-scale structure of the Lya 
forest. A cross-correlation analysis of these data, performed 
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by [Williger et al. (199£ ), revealed a 3.7a excess of l ine-pairs 
at velocity splittings 50 < An < 100 km s _1 . Thus Williger 
et al.lcdmcluded that the Lya forest seems to exhibit struc- 
ture on the scale of ~ 10 h" 1 Mpc in the plane of the sky. 
However there was no excess at smaller velocity splittings 
and no dependence of the signal on angular separation was 
found. 

Finall y, the largest line of sight separa tions were inves- 
tigated by Elowitz, Green, & Impey (1995), who studied a 
group of four QSOs, projected within 2. 8 on the sky. Prob- 
ing scales of ~ 30 h _1 Mpc in the plane of the sky, no sig- 
nificant cross-correlation signal was found out to a velocity 
separation of 10 4 km s _1 . 

From the results outlined above it appears that the Lya 
forest shows significant correlations across lines of sight at 
all but the largest scales probed. In this paper we identify 
more precisely the upper limit of these cor relations. To th is 



Table 1. QSOs analysed. 



end we re-analyse the group of ten QSOs of Williger et al 



the South Galactic Pole region. We employ a novel technique 
which is based on the statistics of the transmitted flux rather 



than t he statis tics of fitted absorption lines (Liske, Webb, 
Carswc '. 



1 1998 ) (hereafter paper I) . The South Galactic Cap 



region has one of the highest known QSO densities in the sky 
and the dataset is the most useful for large-scale structure 
studies of the Lya forest published so far. We find evidence 
for a transition from strong correlations for proper line of 
sight separations]!] < 3 h _1 Mpc to a vanishing correlation 
for line of sight separations > 6 h _1 Mpc. 

The rest of the paper is organised as follows: in Section 
we briefly describe the data. Section H outlines the method 
of analysis and its main advantages over a 'traditional' two- 
point correlation function analysis. Our results are presented 
in Section and discussed in Section 0. 



I 



2 THE DATA 

The data were gathered as part of a larger survey designed 
to rev eal the large-scale c lustering properties of C IV sys- 
tems ( Williger et al. 1996 ). The observations were made on 
the CTIO 4-m telescope using the Argus multifibre spec- 
trograph. The instrumental resolution was ~ 2 A and the 
signal-to-noise ratio per pixel reached up to 40 per 1 A pixel. 
For a complete description of the ob servations and the re - 
duction process we refer the reader to Williger et al. (1996). 
Of the original sample of 25 QSO spectra only 14 cover 
any part of the region between the Lya and Ly/3 emis- 
sion lines. In order to avoid the proximity effect we exclude 
from the analysis those parts of the spectra which lie within 
3000 km s _1 of the Lya emission line. This leaves us with 
11 spectra. We have excluded an additional spectrum from 
the analysis (Q0043-2606) because of its low S/N and small 
wavelength coverage. We are thu s left with the same sample 
as used by |Williger et al. (1999[ ). We list these QSOs, their 
positions and redshifts in Table \u and show their distribu- 
tion in the sky in Fig. [j|. The data cover 2.17 < z < 3.40. 
The angular separations range from 6.1 to 69.2 arcmin or 



t We use (jo = 0.5, A = and H = 100 h km s" 1 Mpc" 
throughout this paper. 
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Figure 1. Distribution of QSOs in the sky. The field is centered 
on a = 00 h 42 m 10 s and S = -26°40' (B1950). Stars mark the 
positions of the QSOs listed in Table hi and emission redshifts are 
indicated. 



1.4 to 16.3 h x proper Mpc and the emission redshifts range 
from 2.36 to 3.44. 



3 ANALYSIS 

The method of analysis used in this work is described in 
detail in paper I. Our technique does not rely on the statis- 
tics of individual absorption lines but rather identifies local 
over- or underdensities of Lya absorption using the statis- 
tics of the transmitted flux. In order to identify large-scale 
features in a normalised spectrum of the Lya forest we con- 
volve the spectrum with a smoothing function. In fact, we 
smooth the spectrum on all possible scales. On the smallest 
possible scale (= 1 pixel) the spectrum remains essentially 
unchanged whereas on the largest possible scale (= number 
of pixels in the spectrum) the whole spectrum is compressed 
into a single numb er akin to 1 — Da. where Da is the flux 



deficit parameter (|Oke fc Korycansky 1982). When plotted 
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in the wavelength-smoothing scale plane, this procedure re- 
sults in the 'transmission triangle' of the spectrum. The base 
of the triangle is the original spectrum itself and the tip is 



th e values of several ot her parameters: 7 was determined 

by Williger et al. (1999J ) for the present data to be 2.5 and 

we take j3 f rom recent high resolution studie s as 1.5 ( Hrj 

the sinf ;le value which results from smoothing the spectrum et al. 1995| Lu et al. 1996, Kim et al. 1997; Kirkman & 



on the largest possible scale. Thus a cluster of absorption Tytler 1997J). As in paper I we determined the width of the 



lines will be enhanced by filtering out the 'high-frequency 
noise' of the individual absorption lines. We use a Gaussian 
as the smoothing function and denote the smoothed spec- 
trum by G(A, a s ). 

We are interested in identifying local fluctuations 
around the mean Lya transmission and in assessing their 



auto-covariance function of a 'perfect' spectrum (i.e. before 
it passes through the instrument), q, from simulations. We 
simulated each spectrum of the dataset 100 times in accor- 
dance with the null-hypothesis (see also paper I), using the 
7 and P values as above and p,b = 30 km s -1 , at = 8 km s 



statistical significance. To this end we have calculated both 
the expected mean, (G), and variance, a G , of the Lya trans- 
mission as functions of wavelength and smoothing scale in 
paper I. The calculations are based on the simple null- 
hypothesis that any Lya forest spectrum can be described by 
a collection of individual absorption lines whose parameters 
are uncorrelated. Thus, in particular, the null-hypothesis 
presumes an unclustered Lya forest. The derivation also 
uses the observationally determined functional form of the 
distribution of the absorption line parameters, for redshift 
z, column density N, and Doppler parameter 6, 



and b cut = 18 km s' 1 (Hu ct al. 1995| ; [Lu et al. 1996 ; 



Kim 



et al. 1997J JKirkman fc Tytler 199^ . The simulated data 



have the same resolution as the real data but a constant 
(conservative) S/N of 20. The simulations are normalised to 
give the same mean effective optical depth as the real data. 
With the values of all the parameters in place, we can 
transform a given transmission triangle into a 'reduced' 
transmission triangle by 



RG(\,a B ) = 



G-(G) 



OG 



(6) 



77(2, N, b) oc (1 + zYN-' 3 exp 



<b-^f 



2al 



(1) 



1997 



for b > b cut flCarswcll ct al. 1984 |Lu ct al. 1996J ; [Kim et al 



I* irkman & Tytler 1997) and the results are given in 



terms of its parameters: 



<G)(A,a s ) 
and 

a G (X,a s ) 



<4(A) + a c 2 _ T (A) 



2i/7T Os/pS 



2 qHX) +1 



(2) 



(3) 



The reduced transmission triangle shows the residual fluc- 
tuations of the Lya transmission around its mean in terms 
of their statistical significance. 

There are several advantages of the method outlined 
above over a more 'traditional' approach involving absorp- 
tion line fitting and construction of the two-point correlation 
function. Most importantly, our technique does not require 
arduous line-fi tting and side-steps all difficu lties arising from 
line-blending ( Fcrnandcz-Soto et al. 199€ ). In paper I we 
demonstrated that for simulated data our technique is signif- 
icantly more sensitive to the presence of large-scale structure 
than a two-point correlation function analysis, especially at 
intermediate resolution. Moreover, numerical simulations of 



where 

a r -r =T L *^ J ^n 



cosmological structure formation ( Cen et al. 1994 ;_ Miralda- 
Escude et al. 1996; Hernquist et al. 1996J; |Zhang, Amiiiios 



-m- 



Norman 1995; Wadsley & Bond 1997; Theuns et al. 199£ 



a n denotes the noise of a spectrum, ctlsf the width of the 
line spread function, ps the pixel size in A, and A Q = 
1215.67 A. q is the intrinsic width of the auto-covariance 
function of a spectrum (approximated as a Gaussian, see also 
paper I). In paper I we tested the analytic expressions above 
against extensive numerical simulations and found excellent 
agreement. Given tha t the Gunn-Peter son optical depth is 
limited to tgp Ss 0.04 ( Webb et al. 1992 ) and given the num- 
ber density of metal lines, we can safely neglect the effects of 
both of these potential contributors to t c s in equation (|2|) . 
A potential source of error in the method used here is 
its sensitivity to the adopted continuum fit. A small error 
in the zeroth or first order of the fit introduces an arbitrary 
offset from the expected mean transmission. However, this 
is easily dealt with by determining the normalisation of the 
mean optical depth, B, for each spectrum individually: 



indicate that at least the low column density forest arises 
in a fluctuating but continuous medium so that a given ab- 
sorption line may not correspond to a well-defined individual 
'cloud'. Thus a method of analysis based on the statistics of 
the transmitted flux itself is most appropriate and intuitive. 
A high signal-to-noise ratio in the QSO spectra is not 
of great importance for the analysis itself (as the noise is 
quickly smoothed over, see equation 0), but it is important 
for reliable continuum fitting. Uncertainties in the place- 
ment of the continuum have already been dealt with above 
(equation H). The second source of uncertainty are the val- 
ues of the parameters of equa tion (|l|) and we will discuss the 
pertaining effects in Sections 4.2 and i.i . 



B = 



&) 



7+1 



ln[G(A 



c 5 v s,max 



)]• 



(5) 



where u s>ma , denotes the largest possible smoothing scale 
and A c is the central wavelength of the spectrum. Thus we 
fix the normalisation for each spectrum at the tip of its 
transmission triangle. 

Before we can evaluate expressions (El) and rt3j) we need 



4 RESULTS 

4.1 Single lines of sight 

In Fig. we show the result of the analysis described above 
for four of the QSOs listed in Table |lj which show the most 
significant features as described below. The vertical lines 
show the Lya positio ns of known metal sys tems which were 
primarily taken from Williger et al. (1996) (their Table 3). 
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Figure 2. Reduced transmission triangles and normalised spectra of four of the QSOs listed in Table hi The smoothing scale on the 
vertical axis is the FWHM of the smoothing Gaussian. Regions of overdense absorption appear yellow and red, underdense regions appear 
blue and purple. The vertical lines indicate the observed Lya wavelengths of known metal absorption systems. Dashed lines mark the 
positions of metal systems that do not have an associated Lya local minimum (see text). 
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A search using NEDm uncovere d only three addi tional sys- 
tems, all towards Q0042-2627 flYork et al. 199l|). Not sur- 



prisingly we see overdense Lya absorption at the positions 
of all metal systems. We take this as an indication that our 
method correctly identifies overdensities. 

Fig. PJ clearly demonstrates the presence of structures 
on scales as large as many hundred km s" 1 . In all there 
are seven features which are significant at the more than 
— 4(7g level. We now briefly describe these quite significant 
detections: 

Q0041— 2707: This spectrum shows three very sig- 
nificant overdensities in Lya. The first, — 6. 2a a at 
(A,FWHM S ) = (3983 A, 4179 km s^ 1 ), is caused by an 
unsaturated cluster of lines which has no associated metal 
system. The second very significant (— 6.0<tg) overdensity in 
this spectrum lies at (4062 A, 498 km s _1 ) and is caused by 
a saturated blend. It seems to be associated with either or 
both of two metal systems detected in C IV and Si IV, Si II 
and C IV respectively. The third significant system is again 
associated with metals, is significant at the — 5.6ctg level, 
and lies at (4373 A, 745 km s _1 ). 

Q0041— 2658: The most significant overdensity in 
the entire sample is a — 6.3o"g overdensity at (3977 A, 
686 km s _1 ), caused by a broad, clearly blended and sat- 
urated feature which is associated with C II, Si IV, C IV 
and Al II absorption. There is also a — 4.3ctg overdensity at 
(4111 A, 278 km s _1 ), again associated with metals (S II 
and C IV). 

Q0042-2627: A -4.5er G overdensity at (4522 A, 
1071 km s _1 ), with no associated metals, caused by a broad 
cluster of lines. 

Q0043-2633: A -5.2cr G overdensity at (4641 A, 
778 km s~ ), with no associated metals, caused by a broad, 
saturated feature. 



4.2 Local minima 

The probability density function (pdf) of G is inherently 
non-Gaussian, such that e.g. the probability of G lying 
within 3<tg of the mean is smaller than 0.9973. In fact, the 
pdf of G is a function of smoothing scale with the Central 
Limit Theorem guaranteeing a Normal distribution at very 
large smoothing scales. (We have attempted to continually 
remind the reader of this point by the persistent use of the 
notation '<tg'-) In light of this difficulty it is necessary to 
estimate reliable significance levels by using simulated data. 
To further demonstrate and quantify the significance of 
the overdensities presented in the previous section we inves- 
tigate the statistics of local minima in the reduced transmis- 
sion triangles by comparing them with the statistics of local 
minima derived from simulated data. We simply define a 'lo- 
cal minimum' (LM) as any pixel in a reduced transmission 
triangle with 

RG(\,cr s ) < -2.0cr G (7) 

and where all surrounding pixels have larger values. If there 



$ The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 



is more than one LM in the same wavelength bin (but at 
different smoothing scales) then we delete the less significant 
one, since we do not want to count the same structure more 
than once. Given the resolution of the data it is evident that 
unclustered absorption lines will produce LM on the smallest 
smoothing scale (= 1.5 A). We may anticipate that the total 
number and distribution of these LM depend sensitively on 
the parameters of the underlying line distribution (fil) . Thus 
we exclude all LM on the smallest smoothing scale from all 
further analyses in order to reduce the model-dependence of 
our results. 

Note that with this definition, almost all metal lines 
have a Lya LM within 300 km s _1 (vertical lines in Fig. PJ). 

As mentioned in Section ra, we have simulated 100 
datasets (= 1000 spectra), essentially by putting down Voigt 
profiles with parameters drawn from distribution (|l|) (see 
paper I for the exact technique). We then applied the same 
procedures as outlined above: first we constructed the re- 
duced transmission triangles for all simulated spectra and 
then we determined the LM. 

In the real data we found 103 LM whereas the simu- 
lated data yielded on average 51 LM with an rms dispersion 
of 7.6. The maximum number of LM found in the 100 sim- 
ulated datasets was 70. Thus there is an excess of the total 
number of LM over the expected number at the > 99 per 
cent confidence level. In Fig. )a we plot 

iV obs (a;) . 



C(x 



A^exp(x) 



(8) 



where N b B (N exp ) is the observed (expected, as derived from 
the simulations) number of LM and x = RG (significance 
level of LM in units of o G ) or x = FWHM S (FWHM of 
smoothing Gaussian). 

From panel (a) we can see that there is a tendency for 
the excess to be more significant for the stronger LM. For 
the leftmost bin we found two LM in the data but none in 
the 100 simulated datasets. Extrapolating from lower sig- 
nificance levels we find iV cxp = 0.013 for this bin, which is 
the value used in panel (a) . A Kolmogorov-Smirnov test in- 
dicates that the distributions N ^ S {RG) and N C ^ P (RG) dis- 
agree at the 95.6 per cent confidence level. The excess is 
strengthened further and the confidence level is increased to 
> 99 per cent by excluding all LM with smoothing scales 
< 250 km s _1 thus showing that the excess is dominated by 
the larger scale LM. In smoothing scale-space we also ob- 
serve an excess which occurs on scales of up to 1200 km s _1 
(panel b). By excluding LM with RG > —4o~g, the excess 
on scales < 250 km s _1 disappears but it persists on larger 
scales thus showing that it is dominated by the more signif- 
icant LM. 

Q0041-2707 and Q0041-2658 exhibit the most signifi- 
cant str uctu res, as can be seen in Fig. Q (and as discussed in 
Section nAh . If these two objects are removed from the sam- 
ple the excess of LM is slightly reduced but remains signifi- 
cant at the 6.1cr level. The effect also persists if we exclude 
non-significant substructure from the analysis by deleting all 
LM which lie within the subtriangle defined by another more 
significant LM on larger smoothing scales. The removal of 
all LM that lie within 300 km s _1 of known metal lines again 
slightly weakens the excess of LM but does not remove the 
effect. 

If we have used incorrect values for the observational 
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Figure 3. (a) Excess of observed number of local minima in reduced transmission triangles over the expected number (= mean of 100 
simulations) as a function of the significance level of the local minima (see equation H). ±1 rms levels as derived from the simulations 
are indicated by the grey shaded areas, (b) Same as (a) as a function of smoothing scale (FWHM of smoothing Gaussian). The insets 
show the same plots but with expanded ^-scales. 



8 J. Liske et al. 



parameters of equation (llj) then we have over- or underesti- 
mated either the absolute value of the transmission fluctua- 
tions or their statistical significance or both. This may affect 
both the normalisation and the shape of the distributions 
N bs(x). The exclusion of LM on the smallest smoothing 
scale from the analysis was designed to minimise the effect 
of changes of the above parameters on iVexp(a;). We have 
checked the success of this strategy by creating six more 
sets of simulations, each consisting of 10 datasets (= 100 
spectra). For each of these sets we varied one parameter: 
P = 1.7, P = 1.3, B (overall normalisation) decreased and 
increased by 20 per cent, fib = 40 km s _1 , and S/N = 10. 
We have found no significant variation of the total number 
of LM in any of these simulations. The shapes of the distri- 
butions of LM as functions of RG and FWHM S also agree 
well with the distributions for the 'standard' case. The only 
exception was the j3 — 1.3 model which produces less sig- 
nificant LM compared to the 'standard' case (although this 
has little effect on ((RG)). In assessing the robustness of 
the result of Fig. H with respect to our choice of parame- 
ter values we thus only need to consider how the individual 
parameters affect JV b a . 

Given equation (J5J) , ctg is virtually independent of 7 for 
all reasonable values. Thus the only effect of changing 7 is 
to change the size of the fluctuations, G — (G) in equation 
(H). Because we normalise (G) at the central wavelength of 
a spectrum (see equation fl) an increase in 7 will decrease 
the magnitude of overdensities at longer wavelengths while 
increasing the magnitude of overdensities at shorter wave- 
lengths and vice versa. Thus there will only be an apprecia- 
ble effect if overdensities tend to lie to one side of a spectrum 
which is not the case. In addition, a wrong 7 value cannot 
explain why the local minima occur almost exclusively on 
certain velocity scales. Nevertheless, we re-determined the 
LM for the real data for 7 = 2.2 and 7 = 2.8 and the effect 
is small. 

By determining the normalisation of the effective opti- 
cal depth for each spectrum individually we avoid over- or 
underestimation of G — (G) caused by small errors in the 
zeroth and first order of the continuum fits of the individ- 
ual spectra. However, if the continua are systematically low 
or high then we will under- or overestimate the mean opti- 
cal depth normalisation, B, and thus under- or overestimate 
og- However, re-analysing the real data with B increased 
by 20 per cent we still find 98 LM, significantly above the 
expected number of 51 ± 7.6. 

The parameters l ib, at, and b cut a ll affect the parame- 
ter q in equation (g). Kim et al. (1997 ) reported a possible 
evolution of fib with redshift over a range covered by the 
present data. If this were correct and if too small a value for 
q were indeed the reason for the observed excess of LM then 
one would expect the excess to be larger at smaller redshifts 
where fib may be higher. Separating the data into a high and 
a low redshift bin we find N i, B (N c ^ p ) = 53(25 ± 5.0) and 
50(26 ± 5.2) for the low and high redshift bins respectively. 
In any case, we have re-analysed the data for fit — 40 km s _1 
(which corresponds to an increase of q of ~ 30 per cent, since 
q is roughly linear in fib) and still found 84 LM. 

Finally we have considered /3 — 1.3 which also has the 
effect of increasing ac- A re-analysis of the data yielded 
79 LM which is still too large to be compatible with the 
simulations. 



From the tests described above we conclude that the re- 
sult presented in Fig. His quite robust: there exist structures 
in the Lya forest on scales of up to ~ 1200 km s _1 at the 
> 99 per cent confidence level. This result constitutes clear 
evidence for the non-uniformity of the Lya forest on large 
scales despite repeated findings that the two-point correla- 
tion function of ab sorption lines shows no signal on scales 
^ 300 km s" 1 (e.g. ICristiani et al. 1997J). 



4.3 Double lines of sight 

So far we have not taken advantage of the fact that the QSOs 
of Table hi are a close group in the plane of the sky. We shall 
now examine possible correlations across lines of sight. One 
of the advantages of the analysis used above is that it is eas- 
ily applied to multiple lines of sight: transmission triangles 
are simply averaged where they overlap. Given the trans- 
verse separations of the QSOs in this group and assuming 
reasonable absorber sizes, two different lines of sight will not 
probe the same absorbers. According to our null-hypothesis 
of an unclustered Lya forest, two different lines of sight are 
therefore uncorrelated. Thus the variance of a mean trans- 
mission triangle (averaged over multiple lines of sight) is 
given by cr G (\, a B )/n, where n is the number of triangles 
used at (\,<j s ). This procedure should enhance structures 
that extend across multiple lines of sight and surpress those 
that do not. 

Here we present our results for the case n — 2. We have 
sorted all pairs of QSOs into one of three groups according to 
the proper transverse separation, d, of the pair (calculated 
at the redshift of the lower redshift QSO): d < 3 h _1 Mpc, 
3 h _1 Mpc < d < 6 h" 1 Mpc and 6 h" 1 Mpc < d. Table 
H lists those pairs within the groups for which the relevant 
spectral regions o verla p. 

As in Section 4.2 we proceeded to construct (, but this 
time using only the overlap regions of the reduced mean 
triangles. Fig. | shows ((RG) and C(FWHM a ) for the three 
cases listed above. The top row shows the results for small 
transverse separations, the middle row for intermediate, and 
the bottom row for large transverse separations. There is 
evidence for a trend: with increasing separation we detect 
fewer and fewer overdensities at large significance levels and 
large smoothing scales relative to the simulations. 

In panel (a) there is an excess of LM which is some- 
what weaker than that seen for single lines of sight (Fig. Ha) . 
However, we detected one LM at —7.5aa- In contrast, the 
100 simulated datasets revealed no LM with RG < — 6<tg- 
(Extrapolating from lower significance levels gives JV cxp = 
0.0033 for this bin.) This very significant overdensity lies at 
Z — 2.272 and is produced by the near-coincidence in red- 
shift space (Av (a 450 km s _1 ) of the two most significant 
single-line of sight overdensities of the entire sample. These 
two overdensities are found in the spectra of Q0041-2658 
and Q0041-2707. The two lines of sight are separated by 
2.31 h _1 Mpc and are the second closest pair in the sample. 
This is a very clear example of a coherent structure traced 
by Lya absorption extending across two lines of sight. 

In all we detect N obs (N exp ) = 41(14±4.7), 62(28±6.9), 
and 65(25 ± 6.4) LM for small, intermediate, and large sep- 
arations respectively. Thus the excess of the total number of 
LM does not seem to vary. However, the shapes of the dis- 
tributions change quite significantly with line of sight sep- 
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Figure 4. Excess of observed number of local minima in reduced mean transmission triangles of QSO pairs over the expected number 
(= mean of 100 simulations) as a function of the significance level (panels a, c, e) and the smoothing scale (panels b, d, f). The value 
of the leftmost bin of panel (a) is 303. ±1 rms levels as derived from the simulations are indicated by the grey shaded areas. The top 
row shows the result for QSO pairs with d < 3 h _1 Mpc, the middle row for 3 h _1 Mpc < d < 6 h _1 Mpc and the bottom row for 
6 h _1 Mpc < d. 
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Table 2. Pairs of QSOs grouped according to their proper transverse separation, d. 



d<3 h" 1 Mpc 


3 h _1 Mpc < d < 6 h" 1 Mpc 


6 h" 1 Mpc < d 


Q0041-2638 Q0042-2639 


Q0041-2638 Q0042-2627 


Q0041-2638 Q0041-2707 


Q0041-2707 Q0041-2658 


Q0041-2638 Q0042-2656 


Q0041-2638 Q0041-2607 


Q0042-2627 Q0042-2639 


Q0041-2638 Q0043-2633 


Q0041-2638 Q0042-2657 


Q0042-2627 Q0043-2633 


Q0041-2707 Q0042-2714 


Q0041-2707 Q0041-2607 


Q0042-2639 Q0043-2633 


Q0041-2707 Q0042-2656 


Q0041-2707 Q0042-2627 


Q0042-2656 Q0042-2657 


Q0041-2707 Q0042-2657 


Q0041-2607 Q0041-2658 




Q0041-2607 Q0042-2627 


Q0041-2607 Q0042-2714 




Q0041-2658 Q0042-2714 


Q0041-2607 Q0042-2657 




Q0041-2658 Q0042-2657 


Q0042-2627 Q0042-2656 




Q0042-2639 Q0042-2656 


Q0042-2627 Q0042-2657 




Q0042-2639 Q0042-2657 






Q0042-2714 Q0042-2657 






Q0042-2656 Q0043-2633 






Q0042-2657 Q0043-2633 




d = 2.45 h" 1 Mpc 


d = 4.59 h" 1 Mpc 


d = 9.67 b- 1 Mpc 



Table 3. Kolmogorov-Smirnov probabilities that the observed 
and expected distributions of RG values and smoothing scales 
agree for the three line of sight separation groups. 



Pks(RG) Pks(FWHM s ) 



d< 3 

3 < d< 6 

6 <d 



7 x 10"° 
0.011 
0.49 



8 x 10"' 
0.017 
0.46 



aration. In Table ra we list KS probabilities that the simu- 
lated and observed distributions of RG values and smooth- 
ing scales agree for the three groups. At small separations, 
the simulated distributions disagree strongly with the ob- 
servations, producing too few LM at large significance levels 
and at large smoothing scales. However at large separations, 
the distributions agree very well. 

If the observed excess of LM in Figs H and M(a) were 
simply due to an incorrect choice of the values for the pa- 
rameters of equation (hi) as discussed in Section 4.2, then it 
is hard to understand why that excess should be so strongly 
reduced for large line of sight separations. The removal of 
the pair Q0041-2707 - Q0041-2658 leaves the result quali- 
tatively unchanged as does the removal of substructure as 
outlined in Section 4.2. The removal of all LM that have a 
metal system within 300 km s _1 in either of the spectra of 
the pair somewhat weakens the excess of panel (a) but does 
not remove the significance of the effect, since we still ob- 
tain P K s(RG) = 9 x 10" 4 and P KS (FWHM S ) = 3 x 10~ 5 at 
small separations but good agreement at large separations 
(cf. Table |). 

A possible explanation for the trend of the decreasing 
excess with line of sight separation could be that the group 
of close QSO pairs is dominated by those QSOs whose spec- 
tra show the most significant overdensities and that these 
QSOs are absent from the group of large-separation pairs. 
However, by inspection of Table we can see that each of the 
groups contains at least eight of the ten QSOs at least once. 
In particular, Q0041-2707 and Q0041-2658 are present in 
all groups and both appear in the group of close QSOs only 
once. 



To further investigate whether the excess of LM seen 
in the double lines of sight (DLOS) is simply due to the 
excess already detected in the single lines of sight (SLOS) 
we attempt to identify the structures in the SLOS that give 
rise to the LM detected in the DLOS: for each LM in a 
DLOS we search for the most significant LM in each of the 
constituent SLOS within 1000 km s _1 (our results do not 
depend sensitively on this value). If the LM in the DLOS 
is due to a structure extending across the two lines of sight 
then one would expect less significant LM in both of the 
constituent SLOS. However, if the LM in the DLOS is due 
to a strong overdensity which does not extend across both 
SLOS then one would expect to find a LM in only one of 
the SLOS. If LM are found in both SLOS, then they may 
be quite dissimilar and one should be of greater significance 
than the LM in the DLOS. For each LM in the DLOS we 
have computed the following quantity: 



RGi 
RGu 



1,2 



where RGi denotes the LM in the two SLOS and RGu that 
of the DLOS. This quantity measures how the SLOS-LM 
compare with the DLOS-LM they produce. In Fig. H wc 
plot as thick solid lines the histograms of this measure for 
the three groups of QSO pairs (top = small separation, bot- 
tom = large separation). We have also computed the same 
histograms for the 100 simulated datasets and display the 
mean and mean ± 1 rms histograms as thin solid lines and 
grey shaded areas respectively. There is a small number of 
cases where no LM can be found in either of the single sight- 
lines and these are excluded from this analysis. 

All panels of Fig. M show a large excess of the observed 
distributions over the simulated ones. This is not overly sur- 
prising since we already know that the total number of LM 
exceeds the expected number at all separations. However, 
at small separations (panel a) the excess is skewed to values 
RG/ RGu — 1 < 0. A KS test shows that the simulated and 
observed distributions disagree at the > 99 per cent con- 
fidence level. This discrepancy disappears for larger line of 
sight separations. 

In addition, we have also counted the number of DLOS- 
LM which are produced by only a single SLOS-LM and 
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Figure 5. Histograms of the indicated quantity for all local min- 
ima detected in the reduced mean transmission triangles of QSO 
pairs (thick lines). The thin lines show the mean histograms com- 
puted from 100 simulations. ±1 rms levels as derived from the 
simulations are indicated by the grey shaded areas. The top row 
shows the result for QSO pairs with d < 3 h _1 Mpc, the middle 
row for 3 h — 1 Mpc < d < 6 h — 1 Mpc and the bottom row for 
6 h — 1 Mpc < d. We also show KS probabilities that the simulated 
and observed distributions agree. 



found that the excess of such cases increases from 2.2a over 
3.4<r to 4.4(j for increasing line of sight separation. 

Again, the removal of certain subclasses of LM (those 
of the Q0041-2707 - Q0041-2658 pair, substructure, those 
associated with metal lines) does not change the results sig- 
nificantly. 

In summary, the results above show that for small line of 
sight separations two SLOS-LM combine to make a DLOS- 
LM of greater significance, whereas for large line of sight 
separations the structures seen in the SLOS are 'diluted'. 
We thus take the anti-correlation of the excess of LM seen 
in the reduced mean transmission triangles of double-lines 
of sight with sight-line separation as strong evidence for the 
existence of structures on scales of up to 3 h _1 proper Mpc 
in the Lya forest. 



4.4 Correlation with metal lines 

Not surprisingly, all the metal systems that were found in 
the present data have associated overdense Lya absorption. 
We have repeatedly remarked in the sections above that the 
removal of LM that are associated with metal lines does 
not change the results qualitatively. Here, we take the oppo- 
site approach: we have repeated the analysis of the previous 
section only for those DLOS-LM that have a metal system 



— within 1000 km 



in either of the constituent SLOS. The 



corresponding figure to Fig. H shows only a marginal trend 
with line of sight separation, with Pks(^G) = 0.091 for 
small separations, which is entirely due to the close pair 
Q0041-2658 and Q0041-2707. The corresponding figure to 
Fig. H and its KS tests also indicate that the DLOS-LM are 
due to chance alignments or single SLOS-LM at all sepa- 
rations, however the numbers are very low and we cannot 
draw any definite conclusions. 

It must be kept in mind that the subdivision into LM 
with and without meta ls is subject to a st rong selection ef- 
fect since th e work by Cowie et al. (1995) and Songaila & 
Cowie (1996) has clearly shown that the census of C IV in 
the present data must be significantly incomplete. Thus the 
distinction made here is more one between high-column den- 
sity vs. low-column density structures rather than between 
metals vs. Lya only. Therefore in the real data we have se- 
lected LM that are due to high-column density structures 
whereas in the simulations we have basically drawn a ran- 
dom sample of LM since in the simulated data the column 
density is unrelated to the presence of metals. 

Nevertheless, it appears that the high-column density 
structures traced by metal lines do not produce overdense 
Lya absorption at the distances probed by this sample, with 
one notable exception, in contrast to the lower column den- 
sity Lya only systems. On the other hand, there is little 
difference between the distributions of smoothing scales for 
those LM that are associated with metal lines and those that 
are not. A KS test gives a probability of 0.3 that the two 
distributions are the same. 



4.5 Triple lines of sight 

QSO triplets can provide important constraints on the shape 
of Lya absorbers. It is possible that the excess of LM seen 
in the DLOS is mainly due to filamentary structures. If this 
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Table 4. Triplets of QSOs grouped according to their proper transverse separation, 



d < 3 h" 1 Mpc 


3 h" 1 Mpc < d < 6 h" 1 Mpc 


6 h _1 Mpc < d 


Q0042-2627 Q0042-2639 Q0043-2633 


Q0041-2638 Q0042-2656 Q0043-2633 
Q0041-2707 Q0042-2714 Q0042-2657 
Q0041-2658 Q0042-2714 Q0042-2657 


Q0041-2638 Q0041-2707 Q0041-2607 
Q0041-2638 Q0041-2607 Q0042-2657 


d= 2.82 h- 1 Mpc 


d = 4.73 h^ 1 Mpc 


d = 8.92 h" 1 Mpc 



were the case then one would expect LM of close triple lines 
of sight (TLOS) to be caused by two, not three SLOS-LM. 
However, if the absorption is sheet-like in nature then one 
would expect to find three similar SLOS-LM per TLOS-LM. 
We have repeated the analysis of Section 4.3 for triple 
lines of sight. Table H lists triplets of QSOs grouped accord- 
ing to their pairwise transverse line of sight separations as 
in Table 0. The triplets of the first two groups form more 
or less equilateral triangles on the sky. However, due to the 
larger extent of the full group in the DEC direction than in 
the RA direction, the triplets of the last group form some- 
what flatter triangles. Unfortunately, our QSO sample is not 
dense enough to provide more than one to three triplets per 
group. 



Assuming that the absorbing structures are expand- 
ing with the Hubble flow, we find that the line of sight 
and transverse correlation scales are roughly comparable 
(1200 km s _1 H _1 (z = 2.81) = 1.6 IT 1 Mpc) with a sugges- 
tion that the absorbers might be flattened in the line of sight 
direction since we still detect somewhat significant correla- 
tions on transverse scales of 4.6 h _1 Mpc. Furthermore, the 
analysis of the only close QSO triplet of the sample showed 
that many coincident absorption features are common to 
only two spectra, perhaps indicating an elongated shape in 
the plane of the sky. However, no firm conclusions can be 
drawn here until more QSO triplets have been analysed. 

Using a comparatively small sample with signi ficantly 
smaller line of sight separations, Crotts (198£) and [Srotts 



Fij . ^ shows the results. Again we see an anti- fe Fang (1998) found a stronger correlation signal for lines 



correlation of the excess of LM with line of sight separation with W > 0.4 A than for weaker lines. In contrast, Williger 



(panels a, c, and e). However, for the smallest separations we 



find that 7 of the 10 identified TLOS-LM are due to only two 
(and not three) SLOS-LM. The vast majority of SLOS-LM 
are of smaller significance than their respective TLOS-LM. 
Therefore it seems that at least in the case of this particular 
QSO triplet the absorption more often extends across only 
two lines of sight than across three. However, the numbers 
are small and thus we cannot draw any definite conclusions 
from this result. 



5 CONCLUSIONS AND DISCUSSION 

We summarise our main results as follows: 

(1) We have analysed the Lya forest spectra of ten 
QSOs at (z) = 2.81 contained within a ~ 1 deg 2 field using 
a new technique based on the statistics of the transmitted 



flux. Cc mparison with two-point correlation function analy- 
ses (Section |4.2[ ), along with the results of paper I, suggests 
that this new method is more sensitive to the presence of 
large-scale structure than the two-point correlation function 
of individually identified absorption lines. 

(2) We find structure on scales of up to 1200 km s _1 
along the line of sight and on scales of up to 17 h _1 Mpc 
(comoving) in the transverse direction. We confirm the ex- 
istence of large-scale stru cture in the Lya forest at the > 99 
per c ent confidence level (Pando & Fang 1996; Williger et al 



199S) 



(3) We find strong correlations across lines of sight with 
proper separation < 3 h _1 Mpc. For intermediate separa- 
tions the correlation is weaker and there is only little evi- 
dence for correlation at line of sight separations > 6 h _1 Mpc 
(Fig. |4|) . We thus present the first evidence for a dependence 
of the correlation strength on line of sight separation, and 
place an upper limit of 6 h _1 Mpc on the transverse corre- 
lation scale at z = 2.81. 



et al. (1999) found from their analysis of the present SGP 
data that the inclusion of weak lines (W > 0.1 A) strength- 
ened the ir correlation signal. Here we can tentatively con- 
firm the Wlliger et al. (1999) result. This is in agreement 
with the results of Gen & Simcoe (1997) who predicted from 
their numerical simulations that high-column density lines 
have smaller correlation lengths than low-column density 
ones. This is already evident from a visual inspection of 
the three-dimensional distribution of the absorbing gas in 
the simulations at different overdensities. Large overdensi- 
ties are confined to relatively small, more or less spherical 
regions whereas small overdensities form extended filaments 
and sheets. 

It is interesting to compare our results with predictions 
from hydrodynamical simulations. The typical length of the 
low-column density filaments in the s imulations is of the or- 
der of 1 h" 1 Mpc in proper units (Miralda-Escude et al 



1996J ; [Zhang et al. 1998J ). |Cen fe Simcoe (1997] ) performed a 



detailed analysis of double-lines of sight in a ACDM simula- 
tion. They concluded that for proper line of sight separations 
> 500 h _1 kpc any coincident absorption is due to unrelated 
and spatially uncorrelated clouds. They pointed out however 
that a significant amount of power is missing from their sim- 
ulation on the scale of the simulated box size (10 h" 1 Mpc). 
Even after correcting for this effect though, the transverse 
correlation scale predicted from these simulations remains 
significantly smaller (by about a factor of 3) than the one 
derived in this work. To investigate this discrepancy it will 
be necessary to perform a detailed comparison by subjecting 
simulated spectra, drawn from a suitable (i.e. large) simula- 
tion box, to the same a nalys is we have presented here . 

Croft et al. (1998a) and [Nusser fe Haehnelt (1999|) have 



outlined schemes for recovering the power spect rum of mass 



fluctua tions from Lya forest spectra. Recently, Croft et al 



(19981) performed the first such measurement on scales of 
2-12 h _1 Mpc from 19 intermediate-resolution QSO spec- 
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Figure 6. Excess of observed number of local minima in reduced mean transmission triangles of QSO triplets over the expected number 
(= mean of 100 simulations) as a function of the significance level (panels a, c, c) and the smoothing scale (panels b, d, f). ±1 rms levels 
as derived from the simulations arc indicated by the grey shaded areas. The top row shows the result for QSO triplets with pairwise 
separations d < 3 h _1 Mpc, the middle row for 3 h _1 Mpc < d < 6 h _1 Mpc and the bottom row for 6 h — 1 Mpc < d. 
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tra at z = 2.5. Although we have concentrated in this work 
only on identifying typical correlation scales, our results con- 
firm the usefulness of intermediate-resolution data for large- 
scale structure studies when analysing the distribution of the 
transmitted flux directly. Since the Lya forest is thought to 
trace the mass distribution more closely than galaxies we 
are likely to gain the most direct measurement of the bias 
between galaxy and mass clustering by comparin g the power 



spectrum of the Lya forest with that of galaxies ( Croft et al 
1998bQ 
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